Paramagnetic-ferromagnetic transition in a double-exchange model 
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We study paramagnetic - ferromagnetic transition due to exchange interaction between classical 
localized magnetic moments and conduction electrons. We formulate the Dynamical Mean Field 
Approximation equations for arbitrary electron dispersion law, concentration and relation between 
exchange coupling and the electron band width. Solving these equations we find explicit formula for 
the transition temperature T^. We present the results of calculations of the Tc for the semi-circular 
electron density of states. 
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I. INTRODUCTION 



The double-exchange (DE) model [1-3] is one of the 
basic ones in the theory of magnetism. Magnetic order- 
ing appears in this model due to exchange coupling be- 
tween the core spins and the conduction electrons. The 
Hamiltonian of the model is 



H 



— n' C„ 



na/3 



(1) 



where c and c' are the electrons annihilation and creation 
operators, S,i is the operator of a core spin, t„_„' is the 
electron hopping, J is the exchange coupling between a 
core spin and n electrons, & is the vector of the Pauli 
matrices, and a, (3 are spin indices. 

The model (the core spins being treated as classical 
vectors) began to be studied several decades ago [1-3]. 
Because of a recent general interest in manganites, it 
reappeared in the focus of attention (see reviews [4-7] 
and references therein). In particular, application of the 
Dynamical Mean Field Approximation (DMFA) [8] sub- 
stantially advanced our understanding of the properties 
of the DE model [9,10]. Most of the papers dealing with 
the model, starting from classical paper by De Gennes [3], 
considered the DE Hamiltonian with infinite exchange 
(and with the addition of the antiferromagnetic superex- 
change, which is crucial for the explanation of magnetic 
properties of manganites). However lately many proper- 
ties of the model for arbitrary strength of the exchange 
were analyzed [11]. 

In this paper we calculate the temperature of a 
paramagnetic-ferromagnetic transition T^. in a double- 
exchange model for arbitrary electron dispersion law, 
concentration and relation between the exchange cou- 
pling and the electron band width by formulating and 
solving the DMFA equations. We treat the core spins as 
classical vectors. 



II. DMFA EQUATIONS AND Tc 

Spins being considered as classical vectors S„ = m„ 
(with the normalization |m|^ = 1), the problem becomes 
a single electron one. The DE Hamiltonian in a single 
electron representation can be presented as 



Hnn' — t„-n' — JlTln ' CrSn 

Let us introduce Green's function 

G{E) = {E-H)-^ 
and local Green's function 

G\oc{E) = (Gnn{E) 



(2) 



(3) 



(4) 



In the last Equation, the averaging is with respect to ran- 
dom configurations of the core spins. In the framework of 
the DMFA approach to the problem (see [8,10] and ref- 
erences therein) the local Green's function is expressed 
through the the local self-energy E by the equation 



Gxoc[E)^gjE-t{E) 



where 



(5) 



(6) 



is the bare (in the absence of the exchange interaction) 
local Green's function. The self-energy satisfies equation 



G\oc{E) 



1 



G^^l{E) + ^{E) + Jvn-a ^ 



(7) 



where (A(m)) = J X{m)P{m), and P(m) is a probabil- 
ity of a given spin orientation (one-site probability). The 
quantities G and S are 2x2 matrices in spin space. 
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The difference between Eq. (7) and similar equation 
appearing in the Coherent Potential Approximation [12], 
is due to the fact that in our problem the disorder is an- 
nealed (and not quenched). That is the probability of 
a given (random) core spin configuration is determined 
by the energy of electron subsystem interacting with this 
configuration. 

To understand the DMFA assumption allowing to find 
the probability P(m) self-consistently with the solution 
of Eq. (7), notice that Eq. (7) reduces the problem of 
electron scattering by many spins, to scattering by a sin- 
gle spin with the effective potential V = — Jm-cr— S. The 
spin is embedded in an effective medium, described by 
the Hamiltonian + S, and, hence, by the local Green's 
function Gioc- Consider the change in the number of elec- 
tron states in such system due to a single spin. We get 
[12,13] 



AD{E, m) 



-Im In det 



1 + ( Jma + e) Gioc] , (8) 



where the argument of both Gioc and S is E + iO. So the 
change in thermodynamic potential is [14,11,15] 



(9) 



An{m)= f{E)AD{E,ni)dE, 



where ,f(E) is the Fermi function. The DMFA approxi- 
mation for the one-site probability P(m) is: 



P(m) DC exp [-pAn{m)] . 



(10) 



Eqs. (7) and (10) are the system of non-linear (inte- 
gral) equations. However, the system is grossly simpli- 
fied in paramagnetic (PM) phase, where the macroscopic 
magnetization M = 0. In this case P(m) = const, the 
averaging in Eq. (7) can be performed explicitly, and we 
obtain 



9{E) 



9 



2^5-^(£) + S(£)±J' 



(11) 



where E = SI, G = gi — go{E — I])i, where 1 is a unity 
matrix. 

In the ferromagnetic (FM) phase near the Curie tem- 
perature, Eqs. (7) and (10) can be linearized [9,11] with 
respect to M. Thus we reduce the DMFA equations to a 
traditional MF equation [16] 



P(m) oc exp (-3/3TcM • m) . 



(12) 



The parameter Tc is formally introduced as a coefficient 
in the expansion of Ai7(m) with respect to M. Non- 
trivial solution of the MF equation 



M 



(13) 



can exist only for T < Tc, hence Tc is the ferromagnetic 
transition temperature. In all cases Tc turns out to be 



much less the chemical potential, so we can consider elec- 
tron gas as degenerate. In this case 



/CX> pfl 
f{E)...dE= / ...dE, 
- OO J — OC 



(14) 



where the chemical potential (in paramagnetic phase) is 
found from the equation 



n ■ 



2 r 

i-oc 



Im g dE, 



(15) 



and n is the number of electrons per site. 

For the Tc, after straightforward, though lengthy alge- 
bra, we obtain 



2J2 



Im 



9 



(Sg'-g)(l-FSg) _ 2Jfg 

g'+g^ 3 



dE, 



(16) 



where S and g are determined by the properties of the 

system in the PM phase; they are the solutions of Eqs. 

(5) and (11); g' = . Eq. (16) is the main 

result of our paper. 

In the case of weak exchange J — > Eq. (16) takes the 
form 



2J2 
~3~ 



dE. 



(17) 



In fact, this equation is the MF approximation [16] for 
the RKKY Hamiltonian [17], to which the original Hamil- 
tonian (1) can be reduced in the case J <^W. 

In the opposite case J ^ oo we have to consider only 
one spin sub-band. (Because of the symmetry of the 
problem with respect to transformation n — > 2 — n, ev- 
erywhere further on we consider only the case n < 1.) 
After we shift E, S and /i by J, we obtain 



r Im 


' g' + 9^' 


J —OO 


[.9' -252 J 



TT 

where g is found from the equation 
9{E) = ^ ^ 



dE, 



2 g-^{E) + ll{E) 



(18) 



(19) 



III. Tc FOR SEMI-CIRCULAR DOS 

Let us apply Eq. (16) to the case of semi-circular (SC) 
bare density of states (DOS) 




(20) 



Then 



2 



No{s)de 
E-e 



W 




■ (21) 



From Eq. (5) follows 

^ = E-W^g/A-g-\ 
and Eq. (11) can be written as 



5 = |E 



2^^E-W^g/4±J- 



(22) 



(23) 



Eq. (16) takes the form 



J2^r2 rn 



I 

J — ( 



Im 



dE. (24) 



In the case of weak exchange J/W -C n^/^, we obtain 
4J2 



(4y2 _ l)x/T^, 



where 



n = no{y) = 1 + - (sin ^ y + y^/l' 

TT V 



(25) 



(26) 



In the case of strong exchange J/W ^ n^/^, Eq. (18) 
takes the form [10,15] 



2 



67r J_ 



Im 



dE 



WV2 
An 



0^-— tan-V3(l-2/2) 



(27) 



where n = no(y)/2. If, in addition, n <C 1, we obtain 



Tc = —=Wn K, .53 Wx. 
AV2 



(28) 



Notice, that virtual crystal approximation (for J oo) 
[3] gives 



Tc = —Wx w .13 Wx, 
15 



(29) 



thus being imprecise by a factor of 4. 

For arbitrary exchange, the integral in Eq. (24) can be 
calculated only numerically. The results of such calcula- 
tion are presented on FIG. 1. 
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FIG. 1. Tc as a function of electron concentration n for dif- 
ferent relative strcngtlis of the exchange: J/W = .25 (dotted 
line), J/W = 1 (dash-dotted line), J/W = 2 (dashed line), 
and J/W = 20 (solid line). 



IV. NONFERROMAGNETISM 

Note that that our main result (equation for the Tc) in- 
dicates it's own limits of validity. In part of the J/W — n 
plane, Eq. (24) gives Tc < 0. Negative value of Tc means, 
that at any temperature, including T = 0, the para- 
magnetic phase is stable with respect to appearance of 
small spontaneous magnetic moment, and strongly sug- 
gests that the ground state in this part of the phase plane 
is nonferromagnetic [18]. The mechanisms leading to the 
destruction of ferromagnetism can be easily understood 
in the extreme cases of weak and strong exchange. 

It is well known, that for J — > the DE Hamiltonian 
can be reduced to to the RKKY Hamiltonian (see review 
[17] and references to the original papers therein) 



Hrkky 



(30) 



Thus weak exchange between the electrons and the core 
spins generates effective long-range and sign oscillating 
effective exchange interaction between the core spins. 

In the opposite case J ^ oo, we can reduce the DE 
Hamiltonian (2) to effective (classic) t — J model (see 
review [19] and references therein). To do that, first di- 
agonalize the exchange part of the DE Hamiltonian by 
choosing local spin quantization axis on each site in the 
direction of m, the latter being determined by polar an- 
gle 9 and azimuthal angle cj). In this representation the 
Hamiltonian is 



(31) 



where 



ann' = cos — cos — 



, . "n ■ "n' i(A, ,-ih„) 

-h sm -— sm — — e 



■ . 6n' On . 0n' -id> , 

OriTi' = sm — cos — e — cos — sm — e ^" . 
2 2 2 2 



(32) 
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The next step is to apply a canonical transformation 

H^H = e^He-^ = H+[S,H] + ^[S[S, H]] + ... (33) 

which excludes all band-to-band transitions. This can be 
achieved if we chose the operator S in the form 



inn' I d- 



^nn' ^nn' 



2-/ V ^n'n ^n'n J 



(34) 



We have 



[S,H 



bnn' 
K'n 





nn' — ^nn' 

n'n 



[S,[S,H-]]^^ = 2j:i^^„l 1, (35) 

where Inn" = l^nn'' Keeping terms up to the sec- 
ond order with respect to t (and only site-diagonal part 
of the second order terms) we obtain 



Hnn' — J ^ ^mi' ~\~ inn' 



^nn' 
Ojfi'n 



+ ^ Inn" (m„ • m„" -l)a''5„ 



(36) 



In the second quantization the Hamiltonian obtained has 
the form (ignoring the constant term) 



^n ^n' 

COS — cos 

2 2 



, . "n . "n' i(th 

+ sm — sm -—e^^" ' 



2 2 - dUw 

^ ^ Inn' * ^n'^Ti^^n; 



(37) 



where S (d) is the operator of creation (annihilation) of a 
spinlcss electron. Thus strong exchange between the elec- 
trons and the core spins generates effective short-range 
antiferromagnetic exchange between the core spins. 

In conclusion, we formulated the Dynamical Mean 
Field Approximation equations for the double-exchange 
model with classical spins for arbitrary electron disper- 
sion law, concentration and relation between the ex- 
change and the electron bandwidth. In the vicinity of the 
paramagnetic-ferromagnetic transition critical tempera- 
ture Tc, these equations were reduced to a MF equation, 



describing a single spin in an effective field, proportional 
to the macroscopic magnetization, and we obtained ex- 
plicit formula for the T^. For the semi-circular electron 
density of states we explicitly calculated the transition 
temperature as a function of the exchange interaction 
and electron concentration. 
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